Estimation of galactic model parameters in high latitudes with SDSS 
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Abstract: We estimated the Galactic model parameters for a set of 36 high-latitude fields included 
in the currently available Data Release 5 (DR 5) of the Sloan Digital Sky Survey {SDSS), to explore 
their possible variation with the Galactic longitude. The thick disc scaleheight moves from ~550 pc at 
120° < I < 150° to ~720 pc at 250° < I < 290°, while the thin disc scaleheight is as large as ~195 pc in 
the anticenter direction and ~15% lower at \l\ < 30°. Finally, the axis ratio (c/a) of the halo changes 
from a mean value of ~0.55 in the two first quadrants of the Galaxy to ~0.70 at 190° < I < 300°. 
For the halo, the reason for the dependence of the model parameters on the Galactic longitude arises 
from the well known asymmetric structure of this component. However, the variation of the model 
parameters of the thin and thick discs with Galactic longitude originates from the gravitational effect 
of the Galactic long bar. Moreover, the excess of stars in quadrant I (quadrant III) over quadrant IV 
(quadrant II) is in agreement with this scenario. 
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1 Introduction 

The traditional star-count analysis of the Galactic struc- 
ture have provided a picture of the basic structural and 
stellar populations of the Galaxy. E xamples and re 
views of these analyses can be found in|Ba hcalll (loEB) 



iGilmore, Wvse fc Kuiikeiil ll 19891) .iMaie^ ki ( 1993 TlRobin. 
20001') and recentlv lChen et all (|200ll ) and lSieeel et aD 



(2002). The l argest of the obse rvational studies prior 
to the SDSS (|York et al.1 12000| '1 are based on photo- 
graphic surveys. The Basle Halo Program has pre- 
sented the l argest systeinatic photometr ic survey of 
the Galaxy (|Beckejll965l : lFenkartlll9"89al lbll3ldh. The 
Basle Halo Program photometry is currently being re- 
calibrated and re-analysed, using an improve d calibra- 
tion o f the RGU photometric syst em (Buscr fc FenkarO 
fi990l : lBuser, Rong fc Karaalilll998l . ll999l ). More recent 
and future studies are being based on charge-coupled 
device (CCD) survey data. 

Our knowledge of the structure of the Galaxy, as 
inferred from star count data with colour information, 
entered now to the next level of precision with the ad- 
vent of new surveys such as SDSS, 2MASS, CADIS, 
BATC, DENIS, UKIDSS/VISTA, CFH/Megacam, and 
Suprime. Researchers have used different methods to 
d etermine the Galacti c model parameters. In Table 1 
of lKaraah eFall l|2004l ) we can find an exhaustive list of 
the different values obtained for the structural param- 
eters of the discs and halo of the Milky Way. One can 
see directly that there is a refinement in the numerical 
values of the model parameters. The local space den- 
sity and the scaleheight of the thick disc can be given 



as an example. The evaluation of the thick disc have 
steadily mov ed towards shorter sca l eheights, from 1.4 5 
to 0.65 kpc (|Gilmore fc Reidlll983l : [Chen et al.ll200ll ), 
and higher local densities (2-10%). In many studies 
the ra nge of values for th e par ameters is large. Fo r ex- 
ample, fCheniFaD l|200ll ) and lSiegel erahl i|2002l ) give 
H(5cli3§aaG4-fe(i,0%, respectively, for the local space den- 
sity for the thick disc. However, one expects the most 
accurate numerical values for these recent works. That 
is, either the range for the parameters should be small 
or a single value with a small error should be given 
for each of them. It seems that workers have not been 
able to choose the most appropriate procedures in this 
topic. 

Large range or different numerical values for a spe- 
cific Galactic model parameter as estimated by differ- 
ent researchers may be due to several reasons: 1) The 
Galactic model parameters are Gal actic latitude / lon- 
ftude depen dent . The two works of lBuser. Rong fc Karaalil 
19981 . Il999l ) confirm this suggestion. Although these 
authors give a mean value for each parameter, there are 
differences between the values of a given parameter for 
different fields. Also, it has been recently shown that 
the Galacti c model parameter s are Galactic longitude 
depe ndent (iBi ir et a.lJ l2006al lbl: ICabrera-Lavers et al.l 
l2007l: lAk et al.ll2007aH 2) The Galactic model param- 
etcrs are absolute magnitude dependent (jKaraali et al.l 
|2004; Bilir et al. 2006c). Hence, any procedure which 
excludes this argument give Galactic model parame- 
ters spread in a large range. 3) Galactic model param- 
eters change with limiting distance of completeness. 
That is, a specific model parameter is not the same for 
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each set of Galac tic model paramete rs estimated for 
different volumes (jKaraali et al.ll2007D . 

The difference between the Galactic model param- 
eters estimated for fields with different Galactic lati- 
tudes and longitudes can be explained by the influence 
of the disc flaring and warping. The disc of our Galaxy 
is far from being radially smooth and uniform. On 
the contrary, its overall shape presents strong asym- 
metries. While the warp bends the Galactic plane up- 
wards in the first and second Galactic longitude quad- 
rants (0° < I < 180°) and downwards in the third and 
fourth quadrants (180° < I < 360°), the flare changes 
the scaleheight as a function of radial distance. 

T his warp is present in all Galactic component s: 
dust (Drimmel & SperEel"200f'; 'Marshall et al."2006'), 
gas (Burton f988; Drimmel & Spcrgol 200f ; Nakanishi fc So 



l200a :ILe 



Blitz 



Hcilcs 2006; Voskcs fc Burtonl2006i r 
and s tars (jLopez-Corredoira et al.,.2002 : ,Momanv et al.l 
|2006| V All these components have the same node posi- 
tion, and their distributions are asymmetric. However, 
the amplitude of the dust warp seems to be less pro- 
nounced than the stellar and gaseous warps, that s hare 
the same approximate amp litude (Lopcz-Corrcdoir a et al.l 
l2002l : lMomanv et al.l[200g ). 

The stellar and gaseous fiarings for t he Milky Way 
are also compatible (jMomanv et al.l2006t ) , showing that 
hz increases with the galactocentric radius for R > 
5 kpc liC cnt. Dame fc Fazio 1991; Drimmel fc Spergell 
200ll:lNarava n fc Jog 200a : IL6p ez-Corrc doira et al.l2002l : 
Momanv et al. 2006'). The behaviour of this flare in 



I2OO6I : lJuric et al.l 120081 ) whereas the second one is re- 
lat ed to the remnants o f some historical merger events 
(jWvse fc Gilmordl2005l 'l. 

In this paper, we derive the structural parameters 
of two discs and halo of the Galaxy from very recent 
SDSS data to observe possible changes in the param- 
eters with the Galactic longitude. We used about 1.27 
million stars in 36 high-latitude flelds which cover the 
whole longitude interval (0° < I < 360°), and we evalu- 
ated their absolute magnitudes by means of the recent 
procedures which provides accurate distance determi- 
nation. In Sect. 2 we describe the SDSS data, as well 
as the density laws, absolute magnitudes, distances, 
and density functions employed in the analysis. Esti- 
mation of the Galactic model parameters and their de- 
dence with the Galactic longitude is given in Sect. 
3. Finally, our main results are discussed and sum- 
marised in Sect. 4 and 5, respectively. 

2 SDSS 

The SDSS is a large, international collaboration project 
set up to survey 10 000 square-degrees of sky in five 
optical passbands and to obtain spectra of one million 
galaxies, 100 000 quasars, and tens of thousands of 
Galactic stars. The data are being taken with a dedi- 
cated 2.5-m telescope located at Apache Point Obser- 
vatory (APO), New Mexico. The telescope has two in- 
struments: a CCD camera with 30 2048x2048 CCDs 
in the focal plane and two 320 fiber double spectro- 
graphs. The imaging data are tied to a network of 
brighter astrometric standards (which would be satu- 
rated in the main imaging data) through a set of 22 
smaller CCDs in the focal plane of the imaging cam- 
era. An 0.5-m telescope at APO will be used to tie the 
imaging data to brighter photometric standards. 

The SDSS obtains images almost simultaneously 
in five broad bands {u, g, r, i and z)Q cen tred at 3551, 
4686, 6166, 7480 and 8932 A, respectively (|Fukugita et al.l 
II9961 I. The imaging data are automatically processed 
through a series of software pipelines which find and 
objects and provide photometric and astro- 



the central discs of spiral galaxies has not been stud- 
ied so well due to inherent difficulties in separating the 
several contributions to the o bserved counts or flux. 
iLopez-Corredoira et al.l |2003), for example, flnd that 
there is a deflcit of stars compared to the predictions of 
a pure exponential law in the inner 4 kpc of the Milky 
Way, which could be explained as being a flare which 
displaces the stars to higher heights above the plane 
as we move to the Galactic centre. 

In this scenario, where on the one hand the mean 
disc {z — 0) can be displaced as much as 2 kpc between 

the location of the maximum and the minir num ampli- mc aaurc ob ect 

tudes of the warp llDrimmel fc Spergel|[2roi^ : lL6pez-Corredoira gr.alT ru 4- j i 1 f u- 4- 

1 1 1 r ' I 1 1"" I — I'l — metric c alibrations to produce a catalogue of obiects 

l2002l : lMomanv et al.ll20()6l ). and on the other the scale- ^ 
height of the stars can show differences up to 50% of 
the value for hzjRp, ) in the ra nge 5 < R < IQ kpc 
jAlvcs 2000; Lopcz-C orredoira e t al. 200j: lMomanv et all 
I2OO6I ) to flt a global Galactic disc model that accounts 
for all these inhomogeneities is, at the very least, tricky. 
Because of this the results in the Galactic model pa- 
rameters might depend on the sample of Galactic co- 
ordinates used, as the combined effect of the warp 
and flare will be different at different directions in the 
Galaxy, hence at different lines of sight. 

There is an additional reason for explaining the 
differences between the numerical values of a given 
Galactic model parameter estimated in different di- 
rection of the Galaxy, mainly at larger galactocen- 
tric distances. This is the observed overdensity re- 
gions with respect to an axisymmetric halo, where 
two competing scenarios have been proposed for its 
explanation: the fl rst one is concerned with the triax- 
iality of the halo (|Newberg fc Yannvll200^ : IXu et all 



with calibrated magnitudes, positions and structure 
information. The photometric pipeline (jLupton et al.l 
I2OOII ') detects the objects, matches the data from the 
five filters, and measures instrumental fluxes, posi- 
tions, and shape parameters (which allows the classi- 
flcation of objects as "point source", -compatible with 
the point spread function-, or "extended"). The pho- 
tometric calibration is accurate to roughly 2% rms in 
the g, r and i bands, and 3% in u and z, as deter- 
mined by the con s tancy of stellar pop ulation colours 
(|lvezic et al. 1 120041 : iBlanton et a"Lll2005l ). while the as- 
trometric calibration pre cision is better t han 0.1 arc- 
sec rms per coordinate (jPier et al.ll2003l ). The Data 
Rel ease 5 (DR5) imaging catalo gue covers 8000 deg^ 
( Ade lman-McCarthv et aDl2007^ with a detection re- 
peatability complete at a 95% level for point sources 
brighter than the limiting apparent magnitudes of 22.0, 



^Magnitudes in this paper are quoted in the ugriz sys- 
tem to differentiate them from the former one, ugriz 
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22.2, 22.2, 21.3 and 20.5 mag for u, g, r, i and z, re- 
spectively. The data are saturated at about 14 mag in 
g, r and i and about 12 mag in u and z. 

2.1 Observational data and reduc- 
tion 

The data used in this work were taken from SDSS 
(DR5) WEB serveiQ for 36 high-latitude fields (60° < 
b < 65°) covering different Galactic longitude intervals 
(0° < I < 360°). SDSS magnitudes u, g, r, i, and z 
were used in a total number of 2 164 680 stars in 36 
fields equal in size (totally 831 deg^). Although the 
fields are equal in size, (23.08 deg'^) their surface den- 
sities (number of stars per square-degree) are not the 
same, following a specific trend with Galactic longi- 
tude (Fig. [1} . This is the first clue for the dependence 
of the Galactic model parameters with the Galactic 
longitudes. Owing to the SDSS observing strategy, 
stars brighter than go = 14'" are saturated, and star 
counts are not be complete for magnitudes fainter than 
go — 22"". 2. Hence, our work is restricted to the mag- 
nitude range 15 < go < 22 for the evaluation of the 
Galactic model parameters. 

The E{B — V) colour excess was evaluated indi- 
vi dually for each subsample sou r ce by using the maps 
of ISchleeel. Finkbeiner fc David (|l998t ) through SDSS 
query server and this was reduced to total absorption 
Av via Eq. 1 



4000 



Av = 3.1E{B - V). 



(1) 



In order to determine total absorpti ons, Am, fo r the 
SDSS bands, Am/Ay data given by l| 19991 ). i.e. 
1.593, 1.199, 0.858, 0.639 and 0.459 for m = u, g,r,i 
and z were used respectively. Thus, the de-reddened 
magnitudes, with subscript 0, are 



Uo — u — 


Au, 
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zo — z ~ 
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(6) 



The total absorptions Am are available in the SDSS 
query server. 

All the colours and magnitudes mentioned here- 
after will be de-reddened ones. Given that the loca- 
tion of the vast majority of our targets are at distances 
larger than 0.4 kpc, it seems appropriate to apply the 
full extinction from the maps. Actually, when we com- 
bine the distance r = 0.4 kpc (distance 0.35 kpc from 
the Galactic plane) w ith the scaleheight of the dust, 
([Marshall et al.ll2006l . H = 125 pc), we find that the 
total extinction is reduced to 6% of the value at the 
Galactic plane, at the nearest distance of stars in our 
work. 

According to IChen et al.l (|200ll ). the distribution 
of stars in the go/{g — r)o colour-magnitude diagram 
(CMD) can be classified as follows: Blue stars in the 
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^http: / / www.sdss.org/ dr5 / access / index.html 



Figure 1: Star counts at 6=-|-62°.5 for the 36 fields 
available in the DR5. 



range magnitude 15 < go < 18 are dominated by thick- 
disc stars with a turn-ofi^ at {g — r)o ~ 0.33 mag, while 
Galactic halo stars become significant for go > 18 mag, 
with a turn-off at {g — r)o ~ 0.2 mag. Red stars, {g — 
r)o > 1.3 mag, are dominated by thin disc stars at all 
apparent magnitudes. The CMD, go/{g — 7')o, in Fig. 
[5]shows the mentioned populations and, spectral types 
and absolutes magnitude of stars of these populations. 

However, the CMDs and the two-colour diagrams 
for all objects (not presented here) indicate that the 
stellar distributions are contaminated b y extragalactic 
objects as claimed bv lChen et al.l (l200lD . The star/extragalactic 

object separation is based on the "stellarity p arame- 

ter" a s returned from the SExtractor routines (jBertin fc Arnoriti 
Il996l '). This parameter has a value between (high 
extended) and 1 (point source). The separation works 
very well to classify a point source with a value greater 
than 0.8. Needless to say, this separation depends 
strongly on seeing and sky bright ness. We also appl ied 
the "locus-projection" method of lJuric et al.l (|2008l ) in 
order to remove hot white dwarfs, low-redshift quasars, 
and white/red dwarf unresolved binaries from our sam- 
ple. Briefly, this procedure consists of rejecting objects 
at distances larger than 0.3 mag from the stellar locus 

(Fig. m- 

2.2 Density laws 

In this work we adopted the density laws of the Basle 
group (jBuser. Rons fc Karaalll998l . ll999h . Disc struc- 
tures are usually parameterized in cylindrical coordi- 
nates by radial and vertical exponentials: 

Dt{R,z)=nt exp{-\z\/h^,i) exp{- {R - Ro) / h^) , 

(7) 

where z = zq + r-sin(b), r is the distance to the ob- 
ject from the Sun, b the Galactic latitude, zq the 
vertical distance of the Sun from the Galactic plane 
(jJuric et al.ll2008l . 24 pc), R the projection of the galac- 
tocentric distance on the Galactic plane. Rp t he solar 
distance from the Galactic centre (jReidlll993 . 8 kpc), 
hz,i and hi are the scaleheight and scalelength, respec- 
tively, and rii is the normalized density at the solar 
radius. The suffix i takes the values 1 and 2 as long as 
the thin and thick discs are considered. As this study 
focuses on the dependence of the scaleheight and solar 
normalization on the Galactic longitude, we flxed their 
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Figure 2: Colour magnitude diagram for the star 
sample. Spectral types and absolute magnitudes 
are indicated in the horizontal axis. Then thin disc 
stars are dominant at the red peak whereas the 
thick disc (14 < 50 < 18) and the halo {go > 18) 
stars are concentrated at the blue peak. 
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(g-r)o 

Figure 3: {g — r)o/{r — i)o two colour diagram for 
the field centred at I = 90°. Isodensity contours 
show the position of stars at distance d < 0.3 mag, 
from the stellar locus, adopted from I J uric et all 
( 2008( l. Black circles represent the dwarf stars with 
spectral types and absolute magnitudes stated in 
two lines below the diagram. 




scalelengths to 2.4 and 3.5 kpc for thin and t hick discs, 
respectively, according to lJuric et akl l|2008t ) . 

The density law for the spheroid component is pa- 
ra meterized in diffe r ent fo rms. The most common is 
the Ide Vaucouleurd (| 19481 ) spheroid used to describe 
the surface brightness profile of elliptical galaxies. This 
law has been deprojected into three dimensions by 
lYound l| 19761 ) as 

Ds{R)^ns e^p[-7.669{R/Ref^^]/{R/R,f^\ (8) 

where R is the (uncorrected) Galactocentric distance 
in spherical coordinates. Re is the effective radius and 
Us is the normalized local density. R has to be cor- 
rected for the axial ratio k = c/a, 

R=[x' + {z/Kf]'^\ (9) 
z = ZQ+rsmb, (10) 

x=[RI + (z/tanfo)^ - 2Ro{z/ta.nb)cosl]^^^, (11) 

with r the distance along the line of sight and, b and 
I the Galactic latitude and longitude respectively, for 
the fiel d under invest igation. The form used by Basle 
group (jFenkartll 1989al lbll3ldh is independent of effective 
radius but is dependent on the distance from the Sun 
to the Galactic centre, 

Ds{R) = Us exp[10.093(l - {R/ Rof'^)]/{R/ Rof'^ ■ (12) 

2.3 Absolute magnitudes, distances, 
and density functions 

Absolute magnitudes were determined by two different 

procedures. For abs olute magnitudes 4 < Mg < 8 we 

used the procedure of lKaraali. Bilir fc Tuncel l|2005. hereafter KBTl) . 
whereas for 8 < Mg < 10 we adopted the proce dure of 
iBilir. KaraaU fc TungeFI l|2005. hereafter BKTt . The 
cited absolute magnitude intervals correspond to the 
spectral type intervals F0-K5 and K5-M0 respectively. 
In the procedure of KBT, the absolute magnitude off- 
set from the Hyades main sequence, AM^ , is given as 
a function of both {g — r)o colour and 5o.43 UV-excess, 
as follows: 

AM^ = C35o.43 + C2(5o.43 + Cl5o.43 + Co, (13) 

where Sqas is the UV-excess of a star relative to a 
Hyades star of the colour-index {g — r)o = 0.43 which 
corresponds to So.e and which is determined by the 
colour equations between UB V and SDSS photometry 
(KBT). The coefficients Ci (i=0, 1, 2, 3) are functions 
of {g — r)o colour and they are adopted from the work 
of KBT, and where AM^ is defined as the difference 
in absolute magnitudes of a program star and a Hyades 
star of the same {g — r)o colour: 

AMf = Af;-Aff, (14) 

The absolute magnitude for a Hyades star can be eval- 
uated from the Hyades sequence, normalized by KBT 
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(Eq. 15 of KBT). Th i s proc edure is the one used in 
the work of lAk et all (|2007bD . and has two main ad- 
vantages: 1) there is no need to separate the stars 
into different populations, and 2) the absolute magni- 
tude of a star is determined from its UV-excess indi- 
vidually which provides more accurate absolute mag- 
nitudes compared with the procedure 'in-situ', where 
a specific CMD is used for all stars of the same pop- 
ulation. When one uses the last two equations (Eqs. 
ll3l and [Tl)) and Eq. 15 of KBT simultaneously, it gets 
the absolute magnitude Mg of a star. 

The procedure of KBT was defined only for the 
colours 0.09 < (g — r)o < 0.93, which corresponds to 
absolute magnitudes 4 < Afg < 8. Hence, for the ab- 
solute magnitudes 8 < Mg < 10 we used the equation 
of BKT, which provides absolute magnitudes for late- 
type dwarfs: 

Mg = 5.791(g - r)o + 1.242(r- - i)o + 1.412. (15) 

Stars with faint absolute magnitudes are very use- 
ful, since they provide estimates for the space densi- 
ties at short distances relative to the Sun which com- 
bine the space densities for absolutely bright stars at 
large distances, and the local densities of Hipparcos 
(jJahreiss fc Wielenlll997l ). Thus, we have a sample of 
stars with absolute magnitudes 4 < Afg < 10 which 
enables us to evaluate space density functions in the 
heliocentric distances interval 0.4 < r < 25 kpc, corre- 
sponding to a range of distances of 0.4 < 2; < 21 kpc 
from the Galactic plane. This interval is large enough 
to estimate a set of Galactic model parameters and 
test their change with the Galactic longitude. The ab- 
solute magnitudes in question and the corresponding 
spectral types (from early F type to early M type) for 
the locus points in the {g — r)o/{r — i)o two colour dia- 
gram is shown in Fig. [3] for the field centred at / = 90° 
as an example. The local space density in the interval 
4 < Mg < 10 is fiat and it attributes a mean value of 
logarithmic space density D* = 7.49. 

In a conical magnitude-limited volume, the dis- 
tance to which intrinsically bright stars are visible is 
larger than the distance to which intrinsically faint 
stars are visible. The effect of this is that brighter 
stars are statistically overrepresented and the derived 
absolute magnitude s are too faint. T his effect, known 
as Malmquist bias (jMalmauisd Il92(il ) . was formalized 
into the general formula: 

Mg = Mo-a-^^^^^, (16) 

where Afg is the assumed absolute magnitude, Mq is 
the absolute magnitude calculated for any star using 
KBT calibration, a is the dispersion of the KBT or 
BKT calibration, and A{g) is the differential counts 
evaluated at the apparent magnitude go of any star. 
The dispersion in absolute magnitude calibration of 
KBT and BKT is around 0™.25, corresponding an er- 
ror about 10% in photometric distance. We divided 
stars into the absolute magnitude intervals (4,5], (5,6], 
(6,7], (7,8], (8,9] and (9,10], and we applied the Malmquist 
bias to stars in each interval separately. This approach 



provides (relative) uniform space densities which is 
the essential Malmquist bias. Thus, the corrections 
applied to the absolute magnitudes are 0.005, 0.003, 

0. 007, 0.008, 0.012 and 0.012 for the absolute mag- 
nitude intervals cited above. The correction of the 
Malmquist bias was applied to the SDSS photometric 
data used in this work. 

The combination of the absolute magnitude Afg 
and the apparent magnitude go of a star gives its dis- 
tance r relative to the Sun, i.e., 

[<7-Mg]o = 51ogr-5. (17) 

iGilmore. Wvse fc Joned (jl995 ''l quote an error of ~0.2 
dex in the derivation of [M/H] from the [/i?!/ photom- 
etry for F/G stars which leads to a random uncertainty 
of 20% in the distance estimation. One expects larger 
distance errors for later spectral type stars. However, 
(it — g) and (g — r) colours are more accurate than the 
{U — B) and [B — V) colours which compensates this 
excess error for K stars. 

Logarithmic space densities D* = log D -|- 10 have 
been evaluated for the combination of three population 
components (thin and thick discs and halo), for each 
field where 7? = Af/AVi,2; AVi,2 = (7r/180)2(A/3)(ri- 
ri); A denotes the size of the field (23.08 deg^); r\ and 
r2 are the lower and upper limiting distances of the vol- 
ume AVi,2; A*' is the number of stars per unit absolute 
magnitude; r* — [{r\ + ?'2)/2]^''^ is the centroid dis- 
tance of the volume AVi,2; and z* = r* sin(6), h being 
the Galactic latitude of the field centre. The limiting 
distances of completeness, zi, can be calculated from 
the following equations: 

b-Mg]o = 51ogn-5, (18) 
zi = n sin(fe), (19) 

where gi is the limiting apparent magnitude (15 and 
22, for the bright and faint stars, respectively. Fig. 
ri is the limiting distance of completeness relative to 
the Sun, and Afg corresponds the absolute magnitude 
defining the interval (Mi, M2] where (Afi, Af2) is (4,5], 
(5,6], (6,7], (7,8], (8,9] and (9,10]. For the limiting dis- 
tance of completeness at short and large distances Afg 
assumes the bright and faint absolute magnitudes for 
each absolute magnitude interval. That is, the limiting 
distance of completeness is defined for each absolute 
magnitude to a limited range of spectral types. 

We present the distribution of (g — r)o colours and 
[M/H] metallicities for stars in the field centred at 
I = 90°, as an example, to show the variation in these 
parameters as a function of apparent go and absolute 
Mg magnitudes. Fig. [5] shows how for the same ap- 
parent magnitude interval the peaks of the histograms 
move to redder colours when one goes to fainter abso- 
lute magnitudes. This result is in agreement with the 
fact that stars closer to the Sun are late-type stars, 

1. e. thin disc stars. Another result is that the peaks of 
the histograms for stars with fainter apparent mag- 
nitudes (18 < go < 22), but with the same abso- 
lute magnitudes as the brighter ones, o ccupy bluer 
colours. This confirms the suggestion of IChen et al.l 
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Figure 5: The {g — r)o colour distribution as a function of apparent and absolute magnitudes, for the field 
centred at Z = 90°. 
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Figure 6: The metallicity distribution as a function of apparent and absolute 
centred at Z = 90°. 
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Figure 4: Apparent magnitude histogram for point 
sources (white area) and for final stars sample 
(shaded area) for the field centred at / = 90°. 



(j200ll ). who demonstrated in their Fig. 6 that appar- 
ently fainter stars {go > 18) are dominated by blue 
stars, i.e. halo stars. 

Additionally, Fig. [6] shows three peaks (at metal- 
rich, intermediate and metal-poor parts) in the metal- 
licity distribution. However, all of them are not con- 
spicuous in the same panel, as they correspond to dif- 
ferent absolute and apparent magnitude intervals. If 
we compare the upper and lower panels for stars with 
the same absolute magnitudes but different apparent 
magnitudes, we notice that the peaks shift to lower 
metallicities when one goes from relatively bright to 
faint apparent magnitude intervals. As different peaks 
mean different components of our Galaxy, i.e. thin and 
thick discs and halo, the result stated a bove confirms 
our previous finding (jKaraali et al.l[2007h that different 
populations are dominant at different absolute magni- 
tudes. Most of the stars with fainter absolute magni- 
tudes, 6 < Alg < 7 for example, are metal rich stars 
(disc stars), whereas the absolutely brightest stars, i.e. 
4 < Mg < 5, have low metallicities, [M/H] < -1 
dex, and he nce belong to the halo component of our 
Galaxy. As ICarnevI (|1979| ) quoted, the UV-excess at 
the red end is small which limits the accuracy of the 
metallicity estimation (see Fig. [16] in the Appendix). 
However, the general trend of metallicity distribution 
in Fig. |6]does not give the indication of such an effect. 

We acknowledge that in this work we have not 
applied any correction for binarity or giant/subgiant 
stars. However, most of the evolved stars are prob- 
ably rejected automatically due to the limiting ap- 
parent magnitude at the bright end, i.e. go = 15. 
We compared the number of giants and dwarfs with 
apparent magnitudes go > 15 for a reliable confir- 
mation of our argument. We adopted a mean abso- 
lute magnitude Mg — 1.5 for giants and we calculated 
their corresponding distance from the Galactic plane, 
z = 4.5 kpc, for go = 15. We used the local space 
density D* = 6.35 and the scaleheight (H = 650 pc) 
for the thick disc dwarfs and we compared the cor- 
responding space density with that of thick disc gi- 




5 10 15 20 25 

z (kpc) 



Figure 7: Observed (symbols) and evaluated (thick 
solid lines) space density functions combined for 
stars of all three population components: thin disc 
(thin solid line), thick disc (dashed line) and halo 
(dotted line), at four different Galactic longitudes. 
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Table 1: Galactic model parameters for 36 SDSS high-latitude (60° < b < 65°) fields, resuhing from the 
fits of the analytical density profiles. The c;olumns indicate: Galactic longitude (Z), scaleheight of thin 
{hz,i) and thick discs (/iz,2), local space density of the thick disc {n2/n{) and the halo (na/ni) relative to 
the thin disc, axial ratio of the halo (c/a) and star number density of the field {N), reduced chi-square 
minimum (x^j„), and the corresponding probability. 



Thin disc Thick disc Halo 



< I > 
(°) 


hz 1 
(pc) 


hz 2 

(pc) 


(%) 




(%) 




c/a 




N 

{stars / deg^) 


X 


■2 

min 


Probability 





177±9 


634±46 


10.96±1.23 


0. 


14±0, 


,01 


0, 


,60±0, 


,02 


3672±61 


0, 


,517 


0, 


,972 


10 


178±6 


613±28 


12.50±0.90 


0. 


,16±0, 


.01 


0, 


.58±0, 


.01 


3686±61 


0, 


.488 


0, 


.981 


20 


167±8 


586±27 


15.14±1.15 


0. 


,17±0, 


.01 


0, 


.57±0, 


.01 


3527±59 


0, 


.438 


0, 


.991 


30 


187±9 


620±44 


11.67±1.33 


0. 


17±0, 


,01 


0, 


,56±0, 


,01 


3391±58 


0, 


,844 


0, 


.677 


40 


180±6 


595±28 


13.46±1.46 


0. 


,19±0, 


.02 


0, 


.54±0, 


.02 


3370±58 


0, 


.475 


0, 


.984 


50 


172±6 


609±26 


13.30±0.89 


0. 


,16±0, 


.01 


0, 


.58±0, 


.02 


3251±57 


0, 


.537 


0, 


.965 


60 


195±6 


641±35 


10.33±0.86 


0. 


13±0, 


,01 


0, 


,62±0, 


,01 


3008±55 


0, 


,795 


0, 


,742 


70 


181±7 


616±38 


11.07±1.01 


0. 


12±G, 


,01 


0, 


,62±0, 


,01 


2733±52 


0, 


,335 


0, 


,999 


80 


178±7 


598±35 


10.86±1.03 


0. 


14±0, 


,01 


0, 


,58±0, 


,01 


2514±50 


0, 


,399 


0, 


,995 


90 


192±7 


647±52 


7.89±1.02 


0. 


,12±0, 


.01 


0, 


.60±0, 


.02 


2252±47 


0, 


.541 


0, 


.963 


100 


189±5 


638±35 


7.94±0.69 


0. 


,12±0, 


.01 


0, 


.58±0, 


.01 


2063±45 


0, 


.310 


0, 


.999 


110 


186±4 


648±35 


7.11±0.62 


0. 


12±0, 


,01 


0, 


,59±0, 


.01 


1932±44 


0, 


.479 


0, 


.983 


120 


196±7 


667±68 


5.45±0.90 


0. 


,15±0, 


.01 


0, 


.54±0, 


.01 


1807±43 


0, 


.855 


0, 


.662 


130 


184±6 


581±36 


8.75±0.96 


0. 


,16±0, 


.01 


0, 


.54±0, 


.01 


1732±42 


0, 


.468 


0, 


.986 


140 


172±4 


550±23 


10.26±0.70 


0. 


16±0, 


,01 


0, 


,54±0, 


,01 


1698±41 


0, 


,413 


0, 


,994 


150 


183±4 


568±26 


8.71±0.69 


0. 


,16±0, 


.01 


0, 


.54±0, 


.01 


1712±41 


0, 


.450 


0, 


.989 


160 


195±8 


603±64 


7.23±1.22 


0. 


16±0, 


.03 


0, 


.53±0, 


.04 


1697±41 


0, 


.543 


0, 


.962 


170 


195±7 


638±67 


6.35±1.05 


0. 


1 7-i-n 


ni 


0, 


,54±0, 


,01 


1719±41 


0, 


,816 


n 

\J ■ 


71 5 


180 


195±4 


658±42 


6.22±0.62 


0. 


,14±0, 


,01 


0, 


,61±0, 


,03 


1766±42 


0, 


,681 


0, 


,869 


190 


198±4 


646±37 


6.46±0.59 


0. 


,12±G, 


,01 


0, 


,65±0, 


,01 


1822±43 


0, 


,630 


0, 


,912 


200 


199±7 


657±61 


6.35±0.91 


0. 


,13±0, 


.01 


0, 


.67±0, 


.03 


1911±44 


0, 


.541 


0, 


.963 


210 


179±5 


598±31 


9.46±0.77 


0. 


,14±0, 


.01 


0, 


.71±0, 


.02 


2083±46 


0, 


.396 


0, 


.996 


220 


185±7 


594±46 


9.16±1.12 


0. 


13±0, 


,01 


0, 


,76±0, 


,02 


2186±47 


0, 


,656 


0, 


,892 


230 


198±4 


654±36 


7.08±0.59 


0. 


,13±0. 


.01 


0, 


.75±0, 


.01 


2314±48 


0, 


.372 


0, 


.997 


240 


174±4 


603±26 


9.18±0.57 


0. 


,13±0, 


.01 


0, 


.71±0, 


.04 


2264±48 


0, 


.495 


0, 


.979 


250 


189±4 


621±27 


8.09±0.56 


0. 


16±0, 


,01 


0, 


,64±0, 


,01 


2310±48 


0, 


,331 


0, 


,999 


260 


188±6 


658±44 


7.53±0.79 


0. 


12±0, 


,01 


0, 


,72±0, 


,02 


2439±49 


0, 


,538 


0, 


,964 


270 


186±9 


684±80 


6.82±1.23 


0. 


,11±0, 


,01 


0, 


,73±0, 


,03 


2545±50 


0, 


,399 


0, 


,995 


280 


196±4 


716±35 


6.00±0.46 


0. 


llibO, 


.01 


0, 


.73±0, 


.02 


2728±52 


0, 


.392 


0, 


.996 


290 


188±9 


704±71 


6.85±1.09 


0. 


,10±0, 


.01 


0, 


.75±0, 


.02 


2967±54 


0, 


.501 


0, 


.977 


300 


190±7 


679±50 


7.60±0.89 


0. 


11±0, 


,01 


0, 


,73±0, 


,01 


3133±56 


0, 


,495 


0, 


,979 


310 


173±3 


649±16 


9.33±0.42 


0. 


,12±0, 


.01 


0, 


.68±0, 


.01 


3258±57 


0, 


.390 


0, 


.996 


320 


186±8 


669±43 


8.32±0.86 


0. 


,12±0, 


.01 


0, 


.68±0, 


.01 


3465±59 


0, 


.444 


0, 


.990 


330 


178±7 


600±34 


10.74±1.04 


0. 


,07±0, 


,01 


0, 


,66±0, 


,01 


3590±60 


0, 


,564 


0, 


,953 


340 


184±5 


619±14 


10.30±0.82 


0. 


,15±0, 


.01 


0, 


.61±0, 


.01 


3554±60 


0, 


.398 


0, 


.995 


350 


200±6 


658±39 


8.73±0.89 


0. 


,17±0, 


.01 


0, 


.58±0, 


.01 


3689±61 


0, 


.732 


0, 


.817 


360 


177±9 


634±46 


10.96±1.23 


0. 


14±0, 


.01 


0, 


.60±0, 


.01 


3672±61 


0, 


.517 


0, 


.972 
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Figure 8: The variation of Xmm with three Galac- 
tic model parameters, (a) with scaleheight of the 
thin disc, (b) with scaleheight of the thick disc, 
and (c) with axial ratio of the halo. The figures in 
the panels (a) and (b) are symmetric and the one 
in the panel (c) is a bit skewed to the right. 




1 2 3 4 5 

z (kpc) 

Figure 9: An alternative procedure for the estima- 
tion of the Galactic model parameters, for the field 
centred at I = 90° as an example. In panel (a), 
the space densities based on the observational data 
for distances larger than z = 5 kpc are compared 
with the analytical density law of the halo. Thus, 
only the local space density and the axial ratio of 
the halo are estimated independent of the Galac- 
tic model parameters of the thin and thick discs. 
In panel (b), the density of the halo estimated in 
panel (a) is omitted, and the density function of 
the remaining data for 2: < 5 kpc is compared with 
the combined analytical density laws of the thin 
and thick discs, which provide local space densi- 
ties and scaleheights for the thin and thick discs. 
Finally, in panel (c) only the density function for 
1.5 < 2 < 5 kpc, where the thick disc is dominant, 
is compared with the analytical density law of the 
thick disc. Thus, we estimated the model param- 
eters of the thick disc individually. The Galactic 
model parameters estimated by this procedure are 
in agreement with the ones estimated by compar- 
ison of the combined observational based density 
function with the combined analytical density laws 
for three populations. This agreement excludes 
any possible degeneracy between halo and discs, 
and between two discs. 



10 



Publications of the Astronomical Society of Australia 



ants for local space density D* = 5.58 and the scale - 
height {H = 585 pc) adopted from lBilir et~all (|2006al ). 
It turned out that the number of thick disc giants 
fainter than go — 15™ is less than 8% of the num- 
ber of dwarfs. The number of halo giants relative to 
the number of halo dwarfs at distances z > 4.5 kpc is 
even less, i.e. 2%. We used the selectio n criteria for 
metal-poor giants of iHelmi et all (|2003l ) to estimate 
the cited number of halo giants. These authors de- 
fine the location of the metal-poor giants by the fol- 
lowing criteria: r < 19", 1.1 < (u - g)o < 2.0, 0.3 < 
{g~r)o < 0.8, -0.1 < Pi < 0.6, \s\ > ms-h0.05, where 
Pi = 0.910(M-c/)o-h0.415(g-r)o-1.28, s = -0.249u-h 
0.794g - 0.555r -|- 0.24 and = 0.002. When we ap- 
ply these criteria to stars with 4 < Mg < 10 we obtain 
460 giants corresponding to 2% of the total number of 
stars (N=19 325) for the field 85° < / < 95°. 

The range of the binary stars' fraction is 25-50% 
depending on the spectral types of stars. If we assume 
a binary fraction of 50% then the inferred scaleheight 
in a photometric paralla x evaluation is app roximately 
80% of the actual value l|Siegel et all 120021 ). However, 
as we compare the scaleheight of stars of a specific pop- 
ulation for fields in different directions of the Galaxy, 
we are interested only in the relative values of scale- 
height but not their actual values. Hence, disregarding 
the binarity does not affect our results. 

3 Galactic model parameters 

3.1 Estimation of the Galactic model 
parameters 

We estimated all the Galactic model parameters (the 
local space densities and scaleheights for the thin and 
thick discs, and the local space density and axial ratio 
for the halo) simultaneously, by fitting the space den- 
sity functions derived from the observations (combined 
for the three population components) to a correspond- 
ing combination of the adopted population-specific an- 
alytical density laws. The absolutely faintest stars in 
this work provide space densities at short distances 
from the Galactic plane, z ~ 0.4 kpc. Hence it is pos- 
sible to have reliable extrapolation between them and 
the space density of Hipparcos in the solar neig hbor- 
hood , D* = 7.49 in logarithmic form (Jahreiss & WielenI 
|1997| V corresponding to the mean of the local space 
densities for stars with 4 < Mg < 10. 

We used the classical method for the estima- 
tion of the Galactic model parameters, a method which 
is made in the studies of the Galactic structure that 
has been used in the determination of the most recent 
numerical values for the Galactic model parameters 
(iPhleps et al.l200 0. 2005; Chen et al.l20q^:ISiegel et all 
I2OO2I : iDu et all 12003 . ,20061 : lJuric et al.ll2008l l. The 
comparison of the logarithmic density functions de- 
rived from the observations and the analytical density 
laws are given in Fig. [T] for four fields with Galactic 
longitudes I = 0°, 90°, 180° and 270°. xL^n shows 
almost a symmetrical distribution, as can be observed 
from Fig. |8] which is given as an example. Hence, the 
errors of the Galactic model parameters could be esti- 



mated by changing a given model para meter until an 
incre ase or decrease by 1 was achieved (jPhleps et al.l 
I2OOOI I. All the Galactic model parameters and their 
errors are given in Table [1] The reduced Xmi„ values 
and the corresponding probabilities are also given in 
Table [T] The Xmin values are low, whereas the prob- 
abilities are rather high confirming the reality of the 
Galactic model parameters. 

We also used a different procedure just to test any 
possible degeneracy in the estimation of the model pa- 
rameters: First, we estimated the local space density 
and the axial ratio for the halo by comparing the log- 
arithmic space density function for z > 5 kpc with the 
analytical density law of the halo (Eq. 12). Then, we 
omitted the space density of the halo, estimated by the 
corresponding density law, and compared the new den- 
sity function for z < 5 kpc with the combined density 
laws of thin and thick discs (Fig. |9]). This procedure 
provides Galactic model parameters for thin and thick 
discs independent of the model parameters of the halo. 
The result of the application of the two different pro- 
cedures shows that the corresponding Galactic model 
parameters for a specific population are either identical 
or they diff'er only by a negligible amount. Hence, we 
may argue that no degeneracy exists in the estimated 
parameters. 

A similar procedure is applied to the thin and thick 
discs. We estimated the solar normalizations and the 
scaleheights of thin and thick discs simultaneously by 
the space density function for 0.4 < z < 5 kpc and 
compared the model parameters of the thick disc with 
the corresponding ones estimated by the space density 
function for 1.5 < z < 5 kpc, where the thick disc 
is dominant. Since no significant difi'erences could be 
observed between the compared parameters, we con- 
cluded that no degeneracy exists between the two discs 
either. 

Additionally, we plotted the relative local space 
density of the thick disc versus its scaleheight, and the 
relative local space density of the halo versus its axial 
ratio, for four fields with Galactic longitudes / = 0°, 
90°, 180° and 270° in Fig. [TU]to test the same prob- 
lem. The contours correspond to the same for 
2a and 3o", where a is the standard deviation. In each 
panel, the cross shows the position for the minimum 
X^, which defines the Galactic model parameters on 
the axes with accuracy. 

3.2 Dependence of the Galactic model 
parameters with the Galactic lon- 
gitude 

Table [T] shows that the scaleheights of the thin and 
thick discs, as well as the relative local space densities 
of the thick disc and halo. Even the axial ratio of 
the halo are not the same for 36 fields. That is, these 
Galactic model parameters change as a function of the 
Galactic longitude. 

Fig. [TT] shows the variation in the scaleheight of 
the thin disc (hz) with the Galactic longitude (/). The 
global distribution of hz has a maximum at / « 190°, 
almost in the direction the anti-Galactic centre, whereas 
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Figure 10: Contours of equal obtained for different values of scaleheight and local density for thick 
disc, and axial ratio and local density for halo for the fields with Galactic longitudes of 0°, 90°, 180°, and 
270°. The cross (x) shows the minimum value in each panel while the contours show the la, 2a, and 
3(7 confidence levels. 



the minimum corresponds to the fields in the Galac- 
tic centre direction. However one can separate the 
bins into several sub-samples with segments of differ- 
ent slopes. Additionally, the segments corresponding 
to the fields with Galactic longitudes less than 150° 
have negative slopes whereas the slopes of the ones 
with longitudes greater than 150° are positive. 

The trend of the scaleheight of the thick disc is dif- 
ferent than the one the thin disc (Fig. I12p . The max- 
imum and minimum of the scaleheight for the global 
distribution are / ~ 290° and I « 140°, respectively. 
The slopes of the segments in this figure are also dif- 
ferent than the ones of Fig. 1111 for the same longitude 
set. The error bars are also larger than the ones for the 
scaleheights of the thin disc. Different trends for the 
relative local space densities of the thick disc (712/711) 
and halo (ns/rii), and the axial ratio of the halo (c/a) 
can be also observed in Figs. [T^ and 1131 For exam- 
ple, the minimum of the global distribution of (712/711) 
lies within 180° < I < 200° whereas the minimum of 
(713 /ni) lies in an interval with larger Galactic longi- 
tudes, 270° < I < 290°. We should note that one can 
also observe segments with different slopes in the dis- 
tributions of the relative local space densities for the 
thick disc and halo, and of the axial ratio for the halo. 

Apart from the variation of the Galactic model 
parameters with the Galactic longitude, one can say 
something about the correlation between the estimated 
parameters. Fig. [T5] shows that the scaleheight of the 
thick disc is an increasing function of the Galactic lon- 
gitude in an interval where the local space density of 




140 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
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Figure 11: Variation in the scaleheight of the thin 
disc with the Galactic longitude (the Galactic lat- 
itudes of the fields he within 60° < 6 < 65°). 



the thick disc is a decreasing function and vice versa. 
The same argument holds for the axis ratio and the lo- 
cal space density of the halo (Fig. I13p . Although there 
is a degeneracy between the scaleheight and the local 
space density of the thick disc, and between the axial 
ratio and the local space density of the halo (Fig. llOp . 
the mentioned correlations are real, since the corre- 
sponding reduced xfrUn are rat her low (TableUl . These 
correlations were also cited in lBuser. Rong fc KaraalU 
(1998, 199!|) where seven fields in different directions of 
the Galaxy were investigated with RGU photometry. 

The range, mean, and standard deviations of the 
Galactic model parameters are given in Table [2] The 
mean of the scaleheight of the thin disc is (at least) 
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Figure 12: Variations in the scaleheight and rela- 
tive space density of the thick disc with the Galac- 
tic longitude. 



Table 2: Equations for the variation of the dif- 
ferent Galactic model parameters as a function of 
Galactic longitude (last column). The range, mean 
and the corresponding standard deviation are also 
given in the second, third and fourth columns, re- 
spectively. 



Parameter 


Range 


Mean 


s 


hz,i (pc) 


167-200 


186 


9 


hz,2 (pc) 


550-716 


631 


37 


(n2/ni) (%) 


5.45-15.14 


9.03 


2.34 


(na/m) (%) 


0.07-0.19 


0.14 


0.02 


(c/a) 


0.53-0.76 


0.63 


0.07 



about 30% less than the one that appea red in the hter- 
ature (lOjha et al. 1999: Bu ser. Rong fc Karaah 199^, 
'l999l: lKaraah et al.ll2004t ). For the thick disc, the scale- 
height and the relative local space density are close 
to the ones that appeared in the literature in recent 
years, also obtained from SDSS data. However, the 
upper limit of the space density is a bit higher than the 
cited one up until now (jChen et al ] [200 l]; ISiegel et al.l 
|2002). The axial ratio and the relative local space 
density of the halo are rather close to the correspond- 
ing ones cited by m any authors ('Robin et al.l Il996l: 
Robin. Revle fc Cre zc 2000; Chen ot al. 200ll: ISiegel et al.1 



20021: iKaraaU et al.ll2Q04l: iBiUr et al.1 12006c|) . 



;l ' ' I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' I ' ' t 



0.8 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I [ 




h I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I ri 

30 60 90 120 150 180 210 240 270 300 330 360 
Galactic longitude (degrees) 

Figure 13: Variations in the axis ratio and rela- 
tive space density of the halo with the Galactic 
longitude. 



4 Discussion 

4.1 Evolution of the concept about 
the Galactic model parameters 

Galactic researchers have been working on the mod- 
elling of our Galaxy for about 25 years. The Galactic 
model parameters for the discs and halo have been re- 
fined since this epoch. There is a concensus about the 
idea that the most refined parameters are those of the 
thick disc. Actually, the scaleheight of the thick disc 
decrea sed from the original value of iGilmore fc ReidI 

(| 19831 1. 1.45 kpc, to the recent one, 0.65 kpc whereas 

the s olar normalization increased from 2 to 6-10% ijSiegel et al.l 
|2002| ). Despite the same density laws, different model 
parameters with large ranges ha ve been cited by dif- 
ferent researchers (see Table 1 of lKaraali et alj [2004). 
However, one expects Galactic model parameters ei- 
ther with small errors or with a short range from the 
recent surveys such as SDSS, SEGUE, 2MASS, DE- 
NIS, and UKIDSS which provide accurate magnitude's 
and colours. It seems that we must approach the prob- 
lem from a physi cal point of view. 

We showed in iKaraali et all (|2004l ') that the Galac- 
tic model parameters are absolute magnitude depen- 
dent. The errors for the model parameters estimated 
for a unit absolute magnitude interval are rather small, 
and the numerical values for a specific Galactic model 
parameter increases or decreases, depending on the 
parameter, with the absolute magnitude. For exam- 
ple, the range of the scaleheight of the thin disc is 
264-334 pc for the absolute magnitude interval 5 < 
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Mg < 13. A second example can be given for the 
thick disc. The range of the relative local space den- 
sity of the thick disc for the absolute magnitude in- 
terval 5 < Mg < 9 is 5.25-9.77%, coincident with the 
classical value given in recent works, without paying 
attention to the dependence of this parameter on ab- 
solute magnitude. In other words, the large range of 
the Galactic model parameters is unavoidable in the 
procedure of star counts in which separation of stars 
into different absolute magnitude intervals is not re- 
garded. 

Different absolute magnitude intervals correspond 
to different spectral types and different populations. 
Stars with the brightest absolute magnitudes (4 < 
Mg < 5), intermediate (5 < Mg < 8), and faintest ab- 
solute magnitudes (8 < Mg < 10) are of spectral types 
F-G, G-K and K-M, respectively (Fig. [3. As claimed 
in our previous papers (Karaali et al. 2004; B ilir et aP 
l2006d ) halo, thick disc and thin disc stars are dominant 
at these spectral type intervals, respectively. Hence, 
different Galactic model parameters estimated for dif- 
ferent absolute magnitude diagrams correspond to dif- 
ferent populations and hence stellar ages. 

4.2 Interpretation of the dependence 
of the Galactic model parame- 
ters with the Galactic longitude 

4.2.1 Scenarios for the asymmetric struc- 
ture of the Galaxy and their confir- 
mation by star counts 

The differences between the numerical values of a given 
Galactic model parameter estimated in different di- 
rections of the Galaxy (mainly at large galactocen- 
tric distances), could be explained by the overdensity 
regions with respect to an axisymmetric halo. Two 
competing scenarios have been proposed for this: the 
first one is c oncerned wit h the triaxiality of the halo 
(iNewberg fc Yannv 2006; IXu et all [2009 : lJuric et all 
l2008l ^ ■ whereas the second one is re lated to the rem- 
nants of some his torical events fcf. | Wvse fc Gilmord 
120051 '). Although iNewberg fc Yannvl l|2005l ) claimed 
that the thick disc was symmetric about the Galac- 
tic lon gitude I = 180°; iParker. Humphrevs fc LarsenI 
(|2003l ) interpretered the excess in the numbers of blue- 
and -intermediate- coloured stars above and below the 
Galactic plane in quadr ant I as the asymmetric struc- 
ture o f the thick disc. iParker. Humphrevs fc LarsenI 
(|2003l ') propose similar scenarios for their observations: 
a) fossil remnant of a merger, b) a triaxial thick disc 
or halo, and c) interaction of the thick disc/inner halo 
stars with the bar in the disc. The limiting magnitude 
in our work is much fainter {go — 22) than the limit- 
ing magnitude (0=18) of the authors just mentioned. 
Hence, we thought we can obtain more reliable results, 
if we use a similar procedure. 

W e adopted the procedure of l|Newberg fc Yannvl 
l2005l l. and we plotted the surface density (number of 
stars per square degree of a field) as a function of 
Galactic longitude, for different apparent g- magni- 
tude intervals, i.e. (15, 16], (16, 17], (17, 18], (18, 
19], (19, 20], (20, 21], and (21, 22] (Fig. [g. An 
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Figure 14: Variation of the number of stars with 
different go apparent magnitudes as a function of 
the Galactic longitude. 



asymmetric structure can be observed not only for the 
faint magnitudes which correspond to the halo stars, 
but also for the intermediate magnitudes that favour 
the thick disc stars. All the functions have a flat min- 
imum in the interval 120° < I < 180°, which is more 
conspicuous at fainter magnitudes. That is, there is 
a deficiency in the number of stars in that longitude 
interval. For example, the surface density in this in- 
terval is about 50% less than the one for a field in 
the Galactic centre direction, for the apparent magni- 
tude interval 21 < po < 22. The less surface density 
for the fields in the anti-centre direction relative to 
the fields in the Galactic centre direction is an effect 
of the disc scalelength. However, this effect should 
be most efficient at the field with longitude / = 180° 
which is not the case in Fig. 1141 The excess of stars 
for the fields I ~ 230° corresponds to the "tail of the 
Sagitt arius tidal stre am" overdensity region cited by 
Ne wberg et all (j2006l l. that was also noted bv lXu et al.l 

pooi T 

Following jParker. Humphrevs fc LarsenI (|2003l ). we 
compared the number of stars for five pairs of fields 
which are symmetric relative to the meridian, i.e. the 
plane perpendicular to the Galactic plane and passing 
through the centre of the Galaxy and the Sun. The 
sample pairs are: (30°, 330°), (60°, 300°), (90°, 270°), 
(120°, 240°), and (150°, 210°). The comparison is 
carried out for four colours, i.e. (<? — r)o < 0.35, 0.35 < 
(g - r)Q < 0.60, 0.60 < (g - r)o < 1.20, and 1.20 < 
{g — r)o. As mentioned in the Introduction, blue stars 
in the range 15 < go < 18 are dominated by thick- 
disc stars with a turn-off {g — r)o ~ 0.33, while halo 
stars become significant for go > 18, with a turn-off 
{g — r)o ~ 0.2. Red stars, {g^r)o ~ 1.3 are dominated 
by thin disc stars at all apparent magnitudes. We shall 
treat the ratio of stars given in Table [3] according to 
their colours. We will use the notation Ni/Nj for the 
ratio of stars, for simplicity, where suffixes i and j 
denote the Galactic longitudes of the fields in question. 

(a) The colour range {g — r)o < 0.35 
Apparently faint stars in this interval are halo stars, 
whereas the bright ones are (thin or thick) disc stars. 
The ratio of the number of stars in quadrant I to the 
number of stars in quadrant IV is about 1 or greater 
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Tabk^ 3: Comparison of the surface densities (number of stars per square degree) for fields in different 
quadrants as a function of colour and apparent magnitude. The sufBx denotes the longitude of the field 
in question. 
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than 1 for bright apparent magnitudes, whereas it de- 
creases when one goes to the fainter magnitudes, and 
it approaches about 0.5 at 21 < go < 22. The inequal- 
ity of number of stars for the halo population is due to 
the triaxiality of this component of the Galaxy. The 
same case holds for the halo stars in quadrant II and 
quadrant III. However, for bright stars, the interpre- 
tation of the ratio of stars seems difficult (but see the 
following sections). 

(b) The colour range 0.35 < {g - r)o < 0.60 

This is a perfect colour range for the interpretation of 
the ratios of the number of stars in the corresponding 
fields, since relatively bright stars belong to the thick 
disc population whereas the faint ones are halo pop- 
ulation objects. The ratios N30/N330 and Ne,o/Nsoo 
are greater than 1 for three apparently bright inter- 
vals, i.e. go: (15, 16], (16, 17] and (17, 18] whereas 
it is less than 1 for go: (19, 20], (20, 21] and (21, 
22]. An excess in number of stars in quadrant I rela- 
tive to quadrant IV can be explained by the existence 
of a bar with its near end at the Galactic longitude 
I Ri 27° (jParker. Humphrevs fc Larsen|[2003l '). Accord- 
ing to this argument one expects a deficiency of stars 
in quadrant II relative to quadrant III. Actually, this 
is the case in our work, i.e. N120/N240 and N150/N210 
are less than 1 for (thick) disc stars. Here also, the in- 
equality of number of stars with apparently faint mag- 
nitudes is due to the triaxiality of the halo. The ap- 
parent magnitude interval 18 < g'o < 19 is a transition 
region between disc and halo populations. 

(c) The colour range 0.60 < (5 - r)o < 1.20 

This is also a perfect colour range for discussing the 
number of stars in the fields in different quadrants. 
The ratios A^ao/A'aso and Neo/Nsoo are greater than 
1 for five apparent magnitude intervals i.e. go: (15, 
16], (16, 17], (17, 18], (18, 19] and (19, 20], indicating 
a fainter limiting apparent magnitude for the thick- 
disc stars in this colour range. As expected, the ra- 
tios N120/N240 and N150/N210 are less than 1 for the 
same apparent magnitude intervals. The combination 
of these results confirm our argument stated in the 
previous section, i.e. the inequality of the number of 
stars brighter than go = 20 in quadrant I and quad- 
rant IV originates from the effect of the disc bar. For 
stars with go > 20, all the ratios of the number of stars 
mentioned above are less than 1. These stars belong 
to the population of a triaxial halo, as cited above. We 
should note that the ratio of the number of faint stars 
(halo stars) with this colour range is larger than the 
corresponding ones in the previous colours. 

(d) The colour range {g — r)o > 1.20 

As cited in the previous sections, the red stars belong 
to the thin disc population. Since the number of stars 
is rather small, we omit the brightest apparent magni- 
tude interval, (15, 16]. For the other six intervals, the 
ratios A^ao/A'aao and Nqo/Nsoo are greater than 1, as 
in the previous sections. On the other hand, the ratios 
Afi2o/A'^240 and N150/N210 are less than 1 or equal to 
1. That is, although the asymmetrical structure of the 
disc can also be explained by the data in this section, 
there is a slight difference between the effects of the 
bar on the stars in quadrant I and quadrant HI. 



4.2.2 Relation between the variation of the 
Galactic model parameters and the 
asymmetric structure of the Galaxy 

The global distributions of the Galactic model parame- 
ters in Figs. [TT]-[T3]show minimum and maximum val- 
ues at different Galactic longitudes, whereas the inves- 
tigation of these distributions with a short scale reveal 
interesting sub-structures, i.e. one can separate the 
plots in a figure into several sub-sets which form seg- 
ments with different slopes. It is difficult to interpre- 
tate the dependence of the Galactic model parameters 
by means of such a distribution. However, this can be 
done by combining the observed distributions in Figs. 
Illl - I13l and the results in Section 4.2.1. The difference 
between the corresponding Galactic model parameters 
of the halo for different Galactic longitudes confirm the 
argument that the halo has a triaxial structure. How- 
ever, for the discs the reason is different. Different sur- 
face densities between the fields in quadrants I and IV; 
and in quadrants HI and II confirm the effect of a disc 
bar. The bar induces a grav itational "wake", traps 
and piles up stars behind it plernauist fc WeinbergJ 
ll992l : lDebattista fc Sellwoodlll998l ). Thus, in response 
to a bar, one expects an excess of stars in quadrant I 
(quadrant HI) over quadrant IV (quadrant II) which 
is the case in our work. Also, the disc flarings produce 
a change in the scaleheight of the thin and thick discs 
as the galactocentric distance v aries. T his effect is 
well explained in[L"6 pcz-Corrcdoir a et al.l l(2002. 2004), 
and lCabrera-Lavers et al.l |2007|), who showed that the 
thick disc presents a flaring in the opposite sense as the 
one for the thin disc, that is, a decrease in as R in- 
creases for the thick disc, whe reas there is an increase 
in hz with R for the thin d isc l|L6pez-Corredoira et al.l 
l2002l : lMomanv et al.l[200^ . 

Fig. 1151 shows the trend in both the thin and thick 
disc scaleheights with the galactocentric distance. In 
the upper panel, the solid line represents t he change in 
the sc aleheight of the thin disc obtained bv | L6pez-Corredoira et al] 
(|2002l l. (eq. 4 in ICabrera-Lavers et al.ll2007l ). with /t, 
(Rq) = 185 pc. It is clear that the results obtained 
are compatible with this flaring, but the scatter is too 
high to be conclusive in this statement. Also, the 
range of galactocentric distances covered by the data 
is too small to infer a proper flt from the data. In the 
lo wer panel, the solid line show s the relationship found 
in ICabrera-Lavers et al. I l|2007l . their eq. 5) by using 
2MASS data, at latitudes < b >= 65°, but now consid- 
ermg hz {Rq) = 700 pc0. Data for > 8 kpc are com- 
patible with the trend previously obtained, although 
there is a decrease in the values for the innermost data 
with respect to the predictions of this law that would 
need a further explanation. Therefore, it seems that 
the observed trend in the scaleheight of the thick disc 
via the 2MASS data is also reproducible with the SDSS 
data. Hence the disc flarings have an important con- 
tribution to the observed changes in the parameters 



•^both scaleheights in the solar neighbourhood are de- 
pendent on the populations under consideration, so a slight 
decrease in those values with respect to that obtained for 
the giant population, dominant in the 2MASS data, is not 
unexpected 
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with galactic longitude, as diflterent longitudes imply 
different ranges of galactocentric distances. 

The two points for each < R > bin shown in 
Fig. [15] correspond to the different scaleheights of the 
two symmetric fields relative to the meridian. This 
confirms the asymmetric structure of the thin (upper 
panel) and thick (lower panel) discs. The different be- 
haviours of the radial variation of the scaleheights of 
thin and thick discs allow specific predictions about 
the radial dependence of the vertical velocity disper- 
sion of the discs, since at a given < R> they are in the 
same gravitational field. The different behaviours are 
real, and, as cited at the end of the first paragraph of 
this section, they are due to the opposite sense of the 
flaring of two discs, i.e. decreases as R increases for 
the thick disc, whereas there is an increase in hz with 
R for the thin disc. The variations of the scaleheights 
of thin and thick discs with Galactic longitude are also 
different (Fig. [TT]and the upper panel of Fig. WJ^i . Ad- 
ditionally, one can see more easily that the scaleheights 
of thin or thick discs for two fields that are symmet- 
ric relative to the meridian are different which, again, 
confirms the asymmetrical structure of the two discs. 

Conclusion: The dependence of the Galactic model 
parameters on the Galactic longitude can be explained 
as a result of the combined effects of the triaxial struc- 
ture of the halo, the gravitational effect of the disc bar, 
and the thin and thick disc flarings, which are more 
dominant in some Galactic coordinates with respect 
to others. 
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Figure 15: Variation of the scaleheight of the thin 
(a) and thick (b) discs with < R >. The sym- 
bols (o) and (•) correspond to the data for the 
fields with Galactic longitudes 0° < / < 180° and 
180° < I < 360°, respectively. The error bars for 
the fields with 180° < I < 360° are marked with a 
different symbol (dashed segment) in order to sep- 
arate them from the error bars of the fields with 
0° < ; < 180°. 
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Figure 16: {{u — g)o, {g — r)o) two-colour diagram 
for stars in the field with longitude I ~ 90". The 
difference in UV-excess between the upper and 
lower envelopes of the diagram is small at the red 
end of the diagram which causes an uncertainty 
in the metallicity estimation. The maximum UV- 
excess, A(w — g) = 0.31 mag, corresponds to the 
colour [g ~ r)o = 0.43 mag which is equivalent to 
(B — V)q = 0.60 mag in the photometry. 



